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Abstract 

Introducing the fermionic i?-operator and solutions of the inverse scattering prob- 
lem for local fermion operators, we derive a multiple integral representation for zero- 
temperature correlation functions of a one-dimensional interacting spinless fermion model. 
Correlation functions particularly considered are the one-particle Green's function and 
the density-density correlation function both for any interaction strength and for arbi- 
trary particle densities. In particular for the free fermion model, our formulae reproduce 
the known exact results. Form factors of local fermion operators are also calculated for a 
finite system. 
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1 Introduction 

The evaluation of correlation functions has been one of the challenging problems in research 
on quantum many-body systems in one-dimension, since most of the intriguing phenomena 
induced by underlying strong quantum fluctuations are theoretically described through cor- 
relation functions. In these systems, it is widely known that the existence of models which 
are solvable by means of Bethe ansatz (see [1,2] for example). Though the exact compu- 
tation of correlation functions, of course, is still tremendously difficult even in such models, 
several analytical methods have been recently developed to derive manageable expressions 
for correlation functions, especially in the spin-1/2 Heisenberg XXZ chain. 

In 1990s, correlation functions of the spin-1/2 XXZ chain at zero temperature and for 
zero magnetic field have been expressed as multiple integral forms derived by the g-vertex 
operator approach [3-5]. An alternative method combining the algebraic Bethe ansatz with 
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solutions to the quantum inverse problem has been provided for the XXZ chain in arbitrary 
magnetic fields [6-8] (see also [9] for a recent review). This method can be generalized to the 
finite-temperature and/or the time-dependent correlation functions [10-14]. 

In general, one-dimensional (ID) quantum spin systems are mapped to ID fermion sys- 
tems through the Jordan- Wigner transformation. For the spin-1/2 XXZ chain, the corre- 
sponding system is a spinless fermion model with the nearest neighbor hopping and interac- 
tion. In the thermodynamic limit, the bulk quantities in the XXZ chain are exactly the same 
as those in the spinless fermion model. Namely, in the ID quantum systems, the difference of 
the statistics between the spin and the fermion does not show up, as long as we concentrate 
on their bulk quantities. 

The situation, however, is radically changed when the quantities accompanying a change 
of the number of particles are considered. For instance, let us consider the (equal-time) one- 
particle Green's function (cic| n+1 ) for the spinless fermion model and the transverse spin-spin 
correlation function (^t a m+i) f° r the XXZ chain, which are intuitively the same. In fact, 
due to the difference of the statistics, or equivalently the nonlocality of the Jordan- Wigner 
transformation, both the correlation functions exhibit completely different behavior. For 
instance, one observes an oscillatory behavior for the one-particle Green's function (referred 
to as the A;F-oscillation, where kp is the Fermi momentum), which is peculiar to the fermion 
systems. In contrast, for the transverse spin-spin correlation function, such an oscillatory 
behavior does not appear 1 . 

As mentioned above, exact expressions for the correlation functions of the XXZ chain have 
already been proposed in the form of multiple integral representations. Unfortunately, once 
the Jordan- Wigner transformation is performed and the XXZ chain is considered instead 
of the spinless fermion model, it is difficult to trace the difference of the statistics in the 
framework of multiple integral representations. Namely, to derive a manageable expression 
of correlation functions for the spinless fermion model, we must directly treat the fermion 
system from the very beginning. 

In this paper, introducing the fermionic i?-operator [15] which acts on the fermion Fock 
space, we directly treat the spinless fermion model without mapping to the XXZ chain. 
Combining the method provided in the XXZ chain [7] with solutions to the inverse scatter- 
ing problem of local fermion operators [16], we derive a multiple integral representing the 
equal-time one-particle Green's function and the density-density correlation function at zero 
temperature both for any interaction strengths and for arbitrary particle densities. Our for- 
mulae reproduce the known results for the free fermion model. In addition to the correlation 
functions, we also compute form factors for local fermion operators, which might be useful 
for systematic evaluations of the spectral functions for the spinless fermion model. 

This paper is organized as follows. In the subsequent section, we introduce the fermionic 
i?-operator and the transfer operator constructed by the i?-operator. Then we briefly review 
the algebraic Bethe ansatz for the spinless fermion model. The scalar product of a Bethe state 
with an arbitrary state is presented in section 3. Combining solutions of the inverse scattering 
problem with the scalar product, we compute form factors of local fermion operators. Multiple 
integral representations for correlation functions are derived in section 4. In section 5, using 
the multiple integral representations, we explicitly evaluate correlation functions for the free 

1 In fact, there could exist an oscillatory behavior in the transverse spin-spin correlation function, which, 
however, can be eliminated by a gauge transformation. 
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fermion model. Section 6 is devoted to a brief conclusion. 



2 Spinless fermion model 



The Hamiltonian of the spinless fermion model on a periodic lattice with M sites is defined 
as 

M 



H = H - n c Q - n ?j > 



ra, 



5 " nj+1 ) " I 



(2.1) 



where cj and Cj are the fermionic creation and annihilation operators at the jth site, respec- 
tively, satisfying the canonical anti-commutation relations 



{cj,c k } = {cj,4} = 0, {c],Cfc} = (5 jfc . 



(2.2) 



Here t and A are real constants characterizing the nature of the ground state, and /i c denotes 
the chemical potential coupling to the density operator rij 
Wigner transformation 



CjCj. By use of the Jordan- 



Cj = exp 



j'-i ' 
i7r^n fe 
k=l . 



cj = exp 



j'-i ■ 
-i7r ^ n fc 
fc=i . 



the spinless fermion model (2.1) can be mapped to the spin-1/2 XXZ chain with an external 
magnetic field h: 



M , ^ \ h M 



(2.4) 



where we have changed the variables i and /x c in (2.1) as t — > J and /x c — > /i. 

The difference between the spinless fermion model and the corresponding XXZ chain lies 
only in the difference of boundary conditions. Therefore, in the thermodynamic limit, the 
quantities without accompanying a change of the number of particles such as the density- 
density correlation function (nin m +i) are exactly the same as those such as a longitudinal 
spin-spin correlation function ((1 — cf)(l — a z m+l ))/A for the XXZ chain. As mentioned 
in the preceding section, however, the quantities changing the number of particles such as 
the one-particle Green's function (cic^ , 1 ) exhibit completely different behavior from those 
such as the transverse spin-spin correlation function (o't a m+i)^ ^ ue ^° ^ ne nonlocality of the 
Jordan- Wigner transformation (2.3). Though the exact expression of {&t a m+i) for the XXZ 
chain has already been presented in [7], it is difficult to convert it to the expression for the 
spinless fermion model, as long as we concentrate on the XXZ chain (2.4). Namely we should 
treat the original spinless fermion model without mapping to the XXZ chain. In this section, 
introducing the fermionic -R-operator [15], we directly consider the spinless fermion model 
(2.1). 
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The fermionic -R-operator is defined by 

#12 (A) = 1 - ni - n 2 + . ^ ^ r (ni + n 2 - 2riin 2 ) + -T7TT-^( c i c 2 + 4 c i)> ( 2 - 5 ) 

sh(A + ??) sh(A + rj) 

which acts on V\ <8> s V 2 - Here Vj is a two-dimensional fermion Fock space whose normalized 
orthogonal basis is given by |0)j and := ct|0)j, where Cj\0)j = and (g) s denotes the super 
tensor product. Note that this fermionic i?-operator satisfies the Yang-Baxter equation [15]: 

#i 2 (Ai - A 2 ) J Ri 3 (A 1 )i? 23 (A 2 ) = #23(A 2 )#i3(Ai) J Ri 2 (Ai - A 2 ). (2.6) 

Identifying one of the two fermionic Fock spaces with the quantum space H m , we define the 
L-operator at the mth site by 

L m (X) = R a m{X — £ m ), (2-7) 

where £ m are inhomogeneous parameters assumed to be arbitrary complex numbers. The 
fermionic monodromy operator is then constructed as a product of the L-operators: 

T(A) = L M (A) • • • L 2 (A)Li(A) = A(A)(1 - n a ) + B{X)c a + c\C{\) + D(\)n a . (2.8) 

Thanks to the Yang-Baxter equation (2.6), the fermionic transfer operator defined by 

T(A) = str a T(A) = a (0|T(A)|0) a - a (l|T(A)|l) a = A(X) - D(X), (2.9) 

constitutes a commuting family: [T(A),T(//)] = 0, where the dual fermion Fock space is 
spanned by o (0| and a (l| with a (l| := a (0|c a and o (0|cI = 0. 

The Hamiltonian of the spinless fermion model Hq (2.1) can be expressed as the logarith- 
mic derivative of the fermionic transfer operator T(A) in the homogeneous limit £ m — > rj/2: 



H = tsh(r l )-^\nT(X) 



A = chr/. (2.10) 



Due to the Yang-Baxter equation (2.6), one immediately sees that the following relation is 
valid for arbitrary spectral parameters Ai and A 2 : 

#i 2 (Ai - A 2 )T 1 (A 1 )T 2 (A 2 ) = T 2 (A 2 )T 1 (A 1 )i?i 2 (Ai - A 2 ). (2.11) 

This leads to the commutation relations among the operators A(X), B(X), C(A) and D(X) 
constructing the monodromy operator (2.8). Let us explicitly write down some of these 
relations, which we will use later: 

B(»)B(\) = B(X)B(fi), C(fj)C(X) = C(X)C{fj), 
A(fj)B(\) = f(X,ti)B(X)A(fi) - g(\fi)B(fi)A(X), 
D(fj)B(\) = -/(/x, \)B(\)D(fj.) + gfa, A)B(/i)D(A), 
C(fi)A(X) = /(/x, A)A(A)CV) - gfa, X)A(^)C(X), 
C(»)D(X) = -/(A, M)£>(A)C(/i) + g(X, »)D(»)C(X), 

C^)B(X) + B{X)C{v) = g(X, v)(D(X)A(fi) - D^)A(X)), (2.12) 
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where 

ff\ \ sh(A-^ + r?) shr/ 

/(A '^ )= ACA-/X) ' 5(A '^ ) = ^(A^)- (2 ' 13) 

Let |fi) be the vacuum state for the system (2.1). Obviously \Q) is an eigenstate of the 
transfer operator 

M 

T(»)\n) = (a(»)-d(v)m), o(M) = l, d(M)=n/ _1 (^^)- ( 2 - 14 ) 

rra=l 

A general iV-particle state can be constructed by a multiple action of -B(A) on as 

AT 

1-0) = n ^(A^-)|^>. 7V = 0,1,...,M. (2.15) 

j'=i 

Utilizing commutation relations in (2.12), one finds that the state (2.15) becomes an eigen- 
state of T(/i), if the complex parameters {Xj} satisfy the Bethe ansatz equation; 

( i = d ^ TT ^ Xj ' Afc ) = TT sh ^ Aj ' ~ ^ TT sh ^ Aj ~ Afc + ^ r?i6^ 

1 J " H /(A*, A,) 11 sh(A, - £ m + V ) 11 sh(A, — Afc — 77) ' 1 • ° j 
Then the corresponding eigenvalue t(/z) (i.e. T(/x)|^») = t(p)IV')) is given by 

JV N 

t(m) = aOx) [J f(\j,n) + (-l) N+1 d(n) [J /(/., A,). (2.17) 
i=i i=i 

Note that the existence of an extra factor (— l)^"*" 1 in (2.16) and (2.17) reflects the fermionic 
structure of the system, which does not appear for the XXZ chain (2.4). The relation (2.10) 
leads to the explicit expression of the energy spectrum per site eo of the spinless fermion 
model: 

1 A tsh 2 r/ (I N\ 

^-M^x j + l)M^-l)~^\2-M)- (2 ' 18) 

In the thermodynamic limit, where M — > 00, N — > 00, and N/M is a constant determined 
by a function of the chemical potential /j, c , the Bethe ansatz equation (2.16) characterizing 
the ground state reduces to the following linear integral equation for the distribution function 
of the Bethe roots: 

-27rip tot (A) + J K{X - M )p tot (/x)d/x = t (A, I) . (2.19) 

In the above, the integral kernel K(X) and the function t(X, n) are respectively defined by 

sh2?7 shr; 
( A) = S h(A + r/)sh(A-r/)' * (A ' M) = sh(A-//)sh(A-// + r/)- (2 ' 20) 

The integration contour C = [— A^ c ,A^ c ] is given by e(±A^ c ) = 0, where e(A) is the dressed 
energy satisfying the following linear integral equation: 

-27ri£(A) + j K(X - fi)e(p)dfi = t sh(r/)t (a, |) + /j c . (2.21) 
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Namely, for — 1 < A < 1 (rj = — i£, £ > 0), the contour C is denned by an interval on the real 
axis. In particular, at fi c — > 0, A Mc — > oo. On the other hand, for A > 1 (77 < 0), the Bethe 
roots are distributed on the imaginary axis. In particular, at n c = 0, A^ c = — 7ri/2. Then the 
energy density (2.18) explicitly reads 

e = i shfa) jf t (A, I) p to t(A)dA - Q - (n,)) = jf t (a, |) e(A)dA - ^, (2.22) 

where (n^) is the particle density given by 

(nj> = /" Ptot(A)dA. (2.23) 



c 



Note that the above expression of the ground state energy (2.22), the distribution function 
(2.19) and the dressed energy (2.21) are exactly the same as those for the XXZ chain [1, 2]. 

For later convenience, let us define the inhomogeneous distribution function p(A, £) [7] as 
the solution of the integral equation: 

-27rip(A, O + J K(\ - rfpfji, Odfi = t(X, 0- (2-24) 

Correspondingly, the inhomogeneous total distribution function /9tot(A, {£}) is defined by 

1 M 

Pto 



1 M 

ot(A,m) = -^/>(A,U), (2.25) 



M 

m=l 



which reduces to ptot(A) (2.19), when one takes the homogeneous limit £ m — > rj/2. 

In particular, at zero chemical potential fi c = (corresponding to the half filling case 
(rij) = 1/2), one has 



P(A,0 



2Ceh^(A+^-e) for I A| < 1, C = ir? 

27rZ^n=-oo ch(n»j) ~ 2n I in=l J 4 (j(A+n_£) jg ) ^1 n^l,(/-e 



(2.26) 



where i} n (\,q) is the Jacobi theta function. 



3 Scalar products and form factors 
3.1 Scalar Products 

Our main aim is to evaluate the correlation function of the spinless fermion model (2.1): 

(Wg) 



where \4> g ) is the ground state constructed by substituting an appropriate solution {A} of the 

I? 



Bethe ansatz equation (2.16) into (2.15), and is the dual state: (Vg I = llfci C(Xj). In 
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this paper we concentrate on the one-particle Green's function (Oj = Cj,cj) and the density- 
density correlation function (Oj = rij). Here and in what follows, we consider the "inhomo- 
geneous" model (2.8) for convenience. The correlation functions of the original Hamiltonian 
(2.1) or (2.10) can be obtained by taking the homogeneous limit £ m — > rj/2. 

To apply the algebraic Bethe ansatz to the computation of (3.1), we must express local 
fermionic operators in terms of elements of the monodromy operators (2.8). In fact, such an 
expression has already been obtained by solving the quantum inverse scattering problem for 
local fermion operators [16]. 

Theorem 3.1 [16] For the inhomogeneous spinless fermion model (2.8) with arbitrary in- 
homogeneous parameters £ m (m = 1, . . . M), the local fermionic operators acting on the jth 
space can be expressed in terms of the elements of the fermionic monodromy operator as 

j— i j j— i j 

c) = X\T{t k )B{!i 3 ) \\T-\i k ), Cj = Hn^C^) flr- 1 ^), 

k=l k=l k=l k=l 

j-1 j M 

n 3 = - n n^Ditj) n T- i a k ), i = n n&)- m 

k=l k=l k=l 

To proceed the evaluation by use of the theorem 3.1, one needs to calculate the action 
of the operators A, B, C and D on an arbitrary state (ip\ = (Q\ Y\f=i C(Xj), where {A} is a 
set of arbitrary complex numbers (not necessary the Bethe roots). After simple but tedious 
computation utilizing (2.12), one obtains 

{n\ J] C(Xj)A(X N+1 ) = £ »W 'Ik ' TT^^I n C ^ 

j=i 6=i Uj=i sh(A 7 - X b ) j=1 

N N+1 yiN _X+ V ) N + 1 

(n\Y[c(\j)D{\ N+1 ) =(-i) N d ^> S+i\^ \ W II 
j=i o=i IL=i sh ( Aa " A i) j=i 

3+ a 3+a, 

N N+1 y-fN , \ I \ 

{ n\ n c^) B (x N+1 ) =(-i) n ^ e ^ r^N+i :~ \ f 

j=l a=l Ilk=i sh ( A « - X k) 

k^a 

N+1 . . Uf=l s H X j ~ X a' + V) N+1 

xy — a(V (n\ TT c(\ 3 ). (3.3) 

^(X N+l -X a ,+n) ]J N +\ sh(A, -A a f = \ ' 3 > 

a'+a j^a' j^a,a' 

Note that the extra sign factors (— 1)^ and (— l) N+l appearing in the last two relations 
reflect the fermionic nature of the system (cf. (3.7) and (3.8) in [7]). Finally the action of 
the operator C(A) is free. 

In the above, one distinguishes the terms containing a(Ajv+i) or d(Xjy+i) and the others, 
which are referred to as the direct and indirect terms, respectively. Due to the existence 
of indirect terms, the explicit form of the scalar product (fi| Ylf = i C(fij) Ylf = i B(Xj)\Q) is 
required to compute the expectation value (3.1), where {Xj}f =1 are solutions of the Bethe 
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ansatz equation and {/J-j}jLi are arbitrary complex numbers. This can be explicitly evaluated 
in the same way as in [1, 17]. The difference between the present system and the XXZ chain 
stems from the difference of the commutation relations (2.12). Here we write down only the 
result. 

Proposition 3.1 The scalar product between a Bethe state and an arbitrary state 

N N 

MMIW) = (niHci^llBix^n) (3.4) 

3=1 3=1 



can be expressed as follows: 



-H U.Zb=i sh (Aa -Vb + V) 



N(N- 

n; v <b sh(A a -A 6 ) S h( W 

- Ha) 



mmi{a}) =saK{aiim) = (-i)^— ^ \ , — ^de^Mi™, 



(3.5) 

where {Xj}^ =l are Bethe roots, {[ij}jLi are arbitrary complex parameters. The N x N matrix 
*({/x}|{A}) is defined by 



*i*({M}|{A» = t(Xj,n k ) - (-l) N+1 d(fi k )t(n k , Xj) f[ ZJ^lfy ^ 

and det at denotes the determinant of an N x N matrix. 

Reflecting the fermion statistics, the extra factor (— 1) N+1 enters into the above expression. 
This factor, however, is eliminated by that in the Bethe ansatz equation (2.16), when one 
takes the limit {fi} — > {A}, which gives the square of the norm of the Bethe vector: 

M{A}|{A}) = (-1)^ sh%) fl Sh( ";"\7W ({A}), (3.7) 

a,b=l SnlAa Ab) 

where 



$jfc({A}) = S jk 



d '^ ] X) 



+ K(Xj-X k ). (3.8) 



31 a=l 

Set {n} = . . . , £ n } U {A n+ i, . . . , Xjy} in (3.5). Then the ratio of the determinants of 
^ and $ can be evaluated in the thermodynamic limit: 

^| = n(^tot(A a ,U}))- 1 det n p(A„e fc ), (3.9) 

a=l 

where p(A,£) and Ptot(A,{£}) are defined as (2.24) and (2.25), respectively. 



3.2 Form factors 

Utilizing the scalar product obtained in the preceding subsection, we derive determinant 
representations of the form factors for a finite system. These formulae might be useful in 
calculating the spectral functions of the present model. 
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The form factors of the local fermion operators are defined by 

7V+1 N 

i^HMIW) = (n\ J] c( H )cll[B(x k )\n), 

j=l k=l 
N N+1 

F+(m\{\}\{ri) = (n\ \{C(X k )c m H B( H )\Sl), 

k=l j=l 
N N 

F^(m|{/i}|{A}) = (Q\ [J C( H )(1 - 2n m ) JJ B{X k )\(l), 

3=1 k=l 

\N+1 



(3.10) 



where {X k } k=1 and {//j} -J^ (or {//j}jLi) are solutions of the Bethe ansatz equation (2.16). 
To evaluate them, first we substitute solutions of the quantum inverse scattering problem for 
local fermion operators (3.2) into the above expressions. For instance, F^(m\{[i}\{\}) can 
be written as 

N m-1 M N 

F N (m\{n}\{\}) = 2(n\ n c( H ) n m)A(u) n n - mmiw)- 

3=1 1=1 l=m+l k=l 

(3.11) 

Then, noting that the states ]J B(\ k )\£l) and (Q\ f\ C(fij) are the Bethe states, we remove the 
products of the transfer matrices by (2.17). Finally applying the determinant representation 
of the scalar product (3.5), and using a simple identity Ylf=i TikLi /C\j'j 6c) = 1 derived 
by the Bethe ansatz equation (2.16), we obtain the determinant representations of the form 
factors. They are summarized into the following proposition. 

Proposition 3.2 The form factors F^(m|{/x}|{A}) and F^(m|{A}|{/x}) for local fermion 
operators S m and c m are respectively given by the following determinant representations: 

F N (m\{fi}\{\}) =(-1) 2 

Vm-liiM) U.k=l sh ( A fc - £™ + V) 

det N+1 H-{m\{^}\{\}) 



n; v <v S h( Mfc - H ) n^ </3 s ma q - a^) 

^(-|{A}|M) =-^P^P^^(-IMI{A}), 



Cii({/4)<^ +1 ({/4) 



where the coefficient </>m({A}) is 



N m 

c(W)=nn/( A ^^ 



(3.12) 



(3.13) 



and the (N + 1) x (N + 1) matrix H (m|{/i}|{A}) is defined by 
sh?? 



sh(^ - A A 



HjN+l ~ ti^jiim)- 



N+1 N+1 

J] sh( Ma - A fc + v ) - (-i) N d(\ k ) Yl sh(^ a -x k - v ) 

a=l 



a=l 



/or fc < TV, 
(3.14) 
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On the other hand, the form factor F^(m\{fi}\{\}) for the operator 1 — 2n m is given by 

F N (m\{fx}\{\}) =(-1) 2 „ ■ 

C-i(l A l) J=1 sh ( A i - U + ??) 

x de%(g(M|{A}})-2P(m|Mi{A})) ^ 

nf< fc [sh(A i - A fe ) sh(^ fc - 



where H({/j,}\{\}) is an N x N matrix defined as 



[J sh(^ a - A fe + r?) - (-lJ^+^CAfc) J] sh (^ - Afe - r?) 



(3.16) 



and P(m|{^}|{A}) is an N x N matrix of rank one, 

N 

P jk (m) = t(nj,U) II sh ( Aa " Xk + ( 3 - 17 ) 



a=l 



Here we have set a(A) = 1. Note that, in the derivation of F^(m|{/x}|{A}), we have used the 
orthogonality of the Bethe vectors, and the formula that the determinant of the sum of an 
arbitrary N x TV matrix A and an N x N matrix of rank one B can be given by 

det(^ + ^) = det^ + ydet^), A® = fa M 1 . (3.18) 

^ J for k = I 

Comparing the form factor for the local spin operator of the third component in the 
XXZ chain [18], one finds that the extra factor (— 1) N+1 appears in (3.16). This factor, how- 
ever, is canceled by that in the Bethe ansatz equation (2.16). Hence, in the thermodynamic 
limit where the distribution of the Bethe roots (2.26) characterizing the ground state is the 
same as that of the XXZ chain, F^(m|{/x}|{A}) (3.15) coincides with the form factor of the 
a z m for the XXZ chain 2 . In contrast to this, for the form factor F^(m\{fi}\{\}) (3.12), the 
factor (— 1) N in (3.14) can not be eliminated by the Bethe ansatz equation (2.16). Namely, 
compared with the form factors of the for the XXZ chain, the extra minus sign remains 
in (3.14). This essential difference caused by the statistics between the fermion and the spin 
induces the different behavior of the correlation functions (cic m+1 ) and (ci~cr~ +1 ). 



4 Multiple integral representations for correlation functions 

In this section, we derive multiple integrals representing the one-particle Green's function 
and the density-density correlation function. Inserting the expression of the local fermion 
operators c\ and Cj (3.2) into (3.1), one obtains 



( 1 m+l) " <« g > ~ • ( } 



2 Note that the overall factor (-l)^- 1 )/ 2 in (3.15) and (3.12) can be eliminated by the normalization of 
the Bethe vector (see (3.7)). 
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Note that {c\c m+l ) = -(ci<4 +1 ). 

Following [7], for the density-density correlation function, we conveniently introduce an 
operator Qi >m as Q\. m = YlT=i n k an( ^ cons ider the expectation value of the generating 
function exp(/3<5i,m) {(3 G C): 

wm = ^irc,(»-^)(6.)rc,(A- fl) -(&)w . (4 . 2) 

Then the density-density correlation can be expressed as 

(nm m+1 ) = ^V 2 m -^(eMPQi,m))\ p=Q , (4.3) 

where £> m denotes the lattice derivative defined by P m /(m) = f(m+l)—f(m) and T>^ n f{rn) = 
/(m + 1) — 2/(m) + /(m — 1). In fact, there exist several manners to express the density- 
density correlation by the generating function [7]. For instance, (nin m+ i) is also written by 
using (exp(/3Qi im )n m+ i) as 



d 

(nin m+ i) = V m - 1 — (e-x.p((3Q 1 . m )n m+1 ) 



f3=o' 



(exp(^C^i im jn m+ i) — . (4.4J 

Yrglrg/ 

To evaluate the correlation functions (4.1)-(4.4), we must calculate the multiple action of 
(yl — e@D)(x) on an arbitrary state (ip\ = (Q\ Yif=i where {A} are arbitrary complex 

numbers. Utilizing the commutation relations among A, D and C (2.12), one finds that the 
action on a general state can be written as 

m 

^\H(A-e^D)(x a ) 

a=l 

V n N-n 

= E E Rn({x + }\{x-}\{X^}\{X-})(n\ n C(x+) [I C(A-), (4.5) 

n=0{A}={A+}U{A-} o=l 6=1 

{:r}={:c+}U{:c~} 
|A+|=|x+|=n 

where |A + |, |x + |, etc denote the number of elements of {A + }, {x + }, etc; p = min(m, N). 
Proposition 4.1 The coefficient R n ({x + }\{x~}\{\ + }\{\~}) in (4.5) is given by 

m—n i- n N—n 

Rn({x + }\{x-}\{\+}\{\-})=S n ({x+}\{\+}\{\-}) H a(x-)II/(x+ X") n /(A 6 -,x-) 

a=l L ft=i ft=i 

N-n 

(4.6) 

6=1 6=1 



-e P d(x-)U{-f(x-,x+)} Y[{-f(x-,\;)} 



where the highest coefficient is expressed as 

h-i sh(x+ — Xt + t?) 

^({^ + }|{A + }|{A-}) = na *~ l \ ffl + det TO M jfc (4.7) 

lla<6 [ S HK " A « ) sh ( X « " X b )J 
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with 



N-n 

M jk =a(X+)t(xt,X+)Hf(\-,X+) 

a=l 

N ~ U n ( shQ + - T + -I- 7l) ) 

+ ^ n ( -/<a;, a-,) n{- ^— } \ ■ («) 

The proof is completely parallel to that in the XXZ chain [7] . The difference of the expression 
stems from the difference of the commutation relation of C and D, which reflects the fermionic 
nature of the present system. 

If (^| is a Bethe vector (or equivalently the eigenstate of the transfer operator (2.9)), one 
further has the following expression which is directly applicable to the computation of the 
generating function (4.2). 

Corollary 4.1 If (tp\ = {tt\Y[jLiC(Xj) is a Bethe vector, namely if {X} are the solutions 
to the Bethe ansatz equation (2.16), the coefficient R n ({£, + }\{£,~}\{X + }\{X~}) of the action 
UZM-^^a) is given by 

^({£ + }|{r }|{A + }|{A-}) = 5 n ({£+}|{A+}) 

n m—n N—n m—n N—n n 

x n n fic^b) n n n n/( A «> A ^ m 

a=l 6=1 a=l 6=1 a=l 6=1 

where n + 

£n({£ + }|{A+}) = — ^,6=1 S h (^a ~ A 6 + V) det Jf (4.10) 

n T a <6[sh(A+-A+)sh(e+-4 + )] 1 
Here the n x n matrix M is defined as 

n 
a=l 

Note that here we have set a(A) = 1 and d(^ a ) = 0, and used the fact that the parameters 
{A} satisfy the Bethe ansatz equation (2.16). 

Now we would like to compute the correlation functions. First let us consider the generat- 
ing function (4.2) of the density-density correlation {n\n m+ \). Noting that \tp g ) is the Bethe 
vector characterizing the ground state, one finds that the numerator in (4.2) can be expressed 
W (V'gin^LiC^ — e/3 ^)(Ca)lV ; g}/n^Li r (6i)> where r(A) is the corresponding eigenvalue of 
the transfer matrix (see (2.17)). Then from Corollary 4.1, Proposition 3.1 and the norm of 
the Bethe vector (3.7), the generating function (4.2) is given by a determinant form. By 
d(^ a ) = 0, it immediately follows that the determinant form is the same as that for the gener- 
ating function of the longitudinal spin-spin correlation function (o'iCr^ l+1 ) [7]. Thus utilizing 
the method proposed in [7], we arrive at the following proposition. 



sh(A+-At + r/) sh(A+-g+ + ?? ) 
sh(A+ - A+ + v ) sh(£+ - At + rj) 



(4.11) 
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Proposition 4.2 In the thermodynamic limit M — > oo ; i/te ground state expectation value 
of the generating function exp(/3Qi im ) is expressed as 

dz 3 f jn \ TT TT sh (^a ~ & + */) sh ( A a ~ &>) 



n=0 v ; ,yl j=l ,yc »— i 



C a=l 6=1 Sh ^ a ~ ^ Sh ^ Aa - & + 



where 



x ^ n ({A}|{z})det n [M jfc ({A}|{z})]det„[p(A„z fc )], (4.12) 



and 



XTie contour C depends on the regime and the magnitude of the chemical potential, while T 
encircles and does not contain any other singularities. In the above expression, the 

homogeneous limit £j = r]/2 can be taken trivially. In this case, T surrounds the point rj/2. 

Note that here to replace the sums over the partition of the set {£} with the set of contour 
integrals around {£}, we have used the following relation: 



n m—n 



e nn^>w}) 

ft}=ft+}u{r}«=i 6=i 

|£+|=n 

n , n m ]la=l Y\b=l s H z a ~ Z b ) 

-/ r ndnn/(^^) n:rtnL V-, + ^ «->'' <« 5 » 

where is assumed to be a symmetric function of n variables {z}, and to be analytic 

in the vicinities of Zj = £,k- As for the sum over the partition of the Bethe roots {A} (2.16) 
characterizing the ground state, we have used the following relation, which is valid in the 
thermodynamic limit: 

£ H{x + }) = ^ U n xflM^{mnw), (4.16) 

{A}={A+}U{A-} • JC j=l 

\\+\=n 

where is assumed to be a symmetric function of n variables {A}, and to be zero when 

any two of its variables are the same. 

The multiple integral (4.12) representing the generating function of the density-density 
correlation function is exactly the same as that for the longitudinal spin-spin correlation 
function [7]. The density-density correlation function {nin m+ \) can be obtained by taking 
the derivative with respect to (3 and m as in (4.3). 
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As for the generating function (exp(/?Qi ;m )n m+ i) in (4.4), the numerator can be written 
as (A\ UT=i( A -^ D )(^)D(U+i)\^g)/Ut=i rita). Applying Corollary 4.1, and then using 
the second relation in (3.3), one can obtain a determinant form of the generating function. By 
the same method as in the derivation of (4.12), the following multiple integral representing 
(exp(/3Qi im )n„ t+ i) can be derived in the thermodynamic limit. 

Proposition 4.3 The ground state expectation value of the operator (exp(/3Qi im )n m+ i) in 
the thermodynamic limit is written as 

<exp Wl ,„)n m+l) = £ ' i n p.. I d»«A n n ^ z "~f^f' {X ;~J b \ 

^ ("0 JrfJi 27n Jc „ = i^sh(z a -£ 6 )sh(A a -£ b + 77) 
A sh(A Q -U+i) A sh(A ra+1 -zg + rj) 

x l\ sh(, a -w) M sHx n+1 -x a + v ) Wnamz}) 

x det n [M jk ({X}\{z})]det n+1 [p(X j ,z k ),...,p(X j ,z n ),p(X j ,C m +i)]- (4.17) 

Here the integration contours C and T are taken as in Proposition 4-2. In this representation, 
the homogeneous limit £j = rj/2 can be trivially taken. 

Note that the above expression agrees with that for (exp(/3(5i, m (l— °"m+i)/2)) [7], as expected. 

Finally we would like to consider the equal-time one particle Green's function {c\c^ m+ 1 ) 
(4.1). Because the operators B and C enter into the numerator in (4.1), the computation 
is more complicated than that for the generating functions of the density-density correlation 
function. Nevertheless, applying Proposition 4.1 and then utilizing the third relation in (3.3), 
one finds the method used above is still applicable for the evaluation of (cicj n+1 ). 

Proposition 4.4 The equal-time one-particle Green's function (c\c m+1 ) in the thermody- 
namic limit can be represented as a multiple integral as 



1 I TT dz i [ An+l x f i x n"=l UT=1 S H Za ~ & + V) 



(ci4.+i) = - 



X 



7=2 M*a - 6) 

Ua=i Y\T=2 s H^a ~ &) ilq^i 1 [sh(A n+ i -Zg + n) sh(A TO+2 - z a - rj)] 

n:=jiT=i sh(A a n:=i[sh(An+i -x a+V ) S h(x n+2 -x a - ^ 

x W n ({\}\{z}) de t ra+ iM ifc det ra+2 [p(A j ,zi),... ,p(Aj,z w+ i),p(Aj,g m+ i)], (4.18) 
sn(A„ + i — A n+ 2 + n) 

where the function W n ({X}\{z}) and the (n + 1) x (n + 1) matrix are respectively defined as 
^nUA}|{z}j - — — — ■ — +1 +1 ■ — -, (4.iyj 

Ug=l Ub=l S H^a -h + rj) Ua=l Ub=l sh (Za ~ Z b + 7?) 

anc? 

f/f-x, \A-Lt(\. 7, \ TT n sh(A a -Aj+q) T-rn+l sh(A.,— z 6 +q) , . < 

M jfe = J ( fc ' j) + H A J' Z *Hla=l sh(A 3 -A a+)? ) 116=1 sh( 26 -A 3 +r,) for J < n ^ ^ ^ 

U(>fc,£i) for j = n + 1 
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(a) — 1 < A < 1 

-K 


C A * 






• • 
• • • . 


• • 

2 ' 








- r, C 





(b) A > 1 




Figure 1: The integration contours C and C for — 1 < A < 1 (a) and for A > 1 (b). The 
contour — CuC encircles the singularities {^j}jLi and does not contain any other singularities. 

The integration contour C for the variables {Aj}™^ 1 is taken as in Proposition 1^.2, while 
the contour C for X n+ 2 should be taken such that the contour —C U C surrounds the points 
X n+2 = zi, . . . , z n+ i, £ m+ i which are simple poles of p(X n +2, Zfc) and p(A„ +2 , £ m +i), and does 
not contain any other singularities. The homogeneous limit £j — > rj/2 can be taken trivially 
in the above expression. Correspondingly the contours T and —CL)C encircle the point n/2. 

In figure 1, typical examples of the integration contour C are depicted both for the regimes 
— 1 < A < 1 and A > 1. Recalling that the Fermi point A Mc at zero chemical potential 
(fj, c = 0) is A^ c = oo for — 1 < A < 1 and A Mc = — 7ri/2 for A > 1 (see section 2), one finds 
that the contributions from the right and left edges of C disappear for —1 < A < 1, and those 
from the upper and the lower edges cancel each other for A > 1. Using this together with 
the relation p{\ + rj,z) = — p(A, z) which is valid only for fi c = 0, one arrives at the following 
representation. 

Corollary 4.2 The equal-time one-particle Green's function for the homogeneous case at 
zero chemical potential can be expressed as 

m-l t „ n+1 , „ n+1 /„v,/_ , r;\ \ m 

fi 1 / sh(A a - f) \ m n^i[sh(A w+ i -zg + r,) sh(A ra+2 - z a )] 
U \sh(X a + D) n:=i [sh(A n+1 -X a + v) sh(A n+2 - A a )] 

x Wn({X}\{z}) d et TO+1 M jfc det ra+2 [p(A J , zi), . . . , p(Xj , z n+ ±) , p(Xj , J)], (4.21) 
sn(,A n+ i — A n+ 2j z 

where Mjk is defined as (4.20) for j < n and Mjk = t(zk, t>) f or j = n + 1. 

Compared with the transverse spin-spin correlation function {o~\0~ m+l ) in [7] 3 , one finds the 
sign in front of the second term of Mjk for j < n changes. 

3 There exists a typo in equation (6.13) in [7] for the transverse spin-spin correlation function (o" 1 f a J ^ +1 ). 
It should be corrected by multiplying by —1. 
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5 Free fermion model 



Using the multiple integral representations derived in the preceding section, we reproduce 
the exact expressions of correlation functions for the free fermion model. 

Let us consider the density-density correlation function (nin m +i) first. As mentioned 
before, the generating functions (exp(/3Qi jm )) (4.12) for (nin m +i) is the same as that for 
the spin-spin correlation functions (o'i<J^ n+1 ) for the XXZ chain. Hence the same method 
described in [19] can be applied to derive the exact expression of {nin m+ \). Here we derive 
it by using (4.17) and (4.4). Set rj = -iri/2. Then M jk (4.14) in (4.17) is factorized as 
Mjk = 2(e^ — l)/sh(2(Aj — £&)), which significantly simplifies the integral representation. 
After taking the derivative of (4.17) with respect to (3 and setting = 0, one finds all the 
terms except for n = 1 vanish. Substituting p(X, z) = i/ (it sh(2(A — z)), and taking the lattice 
derivative, one obtains 

/ \ 1 Ia f sh(A 1 -A 2 )y m (zy- m (A 1 ) 

where 

The integral on T can be evaluated by considering the residues outside the contour T i.e. at 
the points z = Ai and A 2 . It immediately follows that 

1 /-, 2x l-y m (A 2 )y- m (A 1 ) 
{niUm+l) = ^J c d ch(2A 1 )ch(2A 2 ) • (5 - 3) 

Changing the variables ch(2A) = — 1/cosA; (k € [kp,2ir — kp]), where kp € [7r/2,7r] is the 
Fermi momentum determined by cos kp = — 1/ ch A Mc (see section 2) , and identifying dA = 
— dfc/(2cosfc), one arrives at 

i />2n — k-p />2tt — fep 

(nm m+1 ) = dk ± dk 2 (l - e- im ^- k ^) 
4vr 2 J k¥ J k¥ 

. , 2 1 - cos(2/c F m) 7r , ,_ A . 

= ^ ~ 2 J m i ^€[-,vr], (5.4) 

where (rij) = 1 — kp/ir is the particle density given by (2.23). This expression agrees with 
the well-known result. At zero chemical potential /x c = 0, where the Fermi momentum 
corresponds to kp = tt/2, one finds 

1 i _ (-l) m 

{ni n m+1 ) = -- 2?r2m2 ■ (5.5) 

Next let us discuss the one-particle Green's function (cic m+1 ) (4.18). Since Mjk = for 
j < n, only the term corresponding to n = survives in (4.18). Namely 

i i 

(5.6) 



Lsh(2(A! - *)) sh(2(A 2 + f i)) sh(2(A 2 - z)) sh(2(Ai + f i))J " 
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To evaluate the above integral explicitly, we shift the A2-contour C to C. Since the integrand 
is antisymmetric with respect to the exchange of Ai and A2, the resulting integral vanishes: 
we only have to take into account the poles surrounded by the A2-contour — C U C i.e. A2 = 
— 7ri/4, z. After evaluating the integral with respect to z by the method used in the calculation 
for the density-density correlation function, one has 

(ClCm+l) = ^Jc dA S h(A + f i)sh(A- f i)- (5 - ?) 
Changing the variables as explained before, one finally obtains 

This also coincides with the well-known result. In particular at zero chemical potential fi c = 
(hp = 7r/2), one sees 

w4 + i> = -(-1 r-^r- ( 5 - 9 ) 

6 Summary and discussion 

In this paper, correlation functions for the spinless fermion model have been discussed by 
using the fermionic i?-operator and the algebraic Bethe ansatz. Applying solutions of the 
inverse scattering problem, we derived determinant representations for form factors of local 
fermion operators. In addition, multiple integrals representing the density-density correlation 
function and the equal-time one-particle Green's function were obtained both for arbitrary 
interaction strengths and particle densities. In particular for the free fermion model, these 
formulae reduce to the known exact results. 

Before closing this paper, we would like to remark possible generalizations of our formulae. 
A generalization to the finite temperature case is of importance. In fact, the quantum transfer 
matrix, which is a powerful tool for study of finite-temperature properties for strongly cor- 
related systems in one-dimension, has already been provided for the spinless fermion model 
[20]. Because the algebraic relations among the elements of monodromy operators in the 
quantum transfer matrix approach are completely the same as (2.12), the relations derived 
in this paper are applicable without essential changes. 

It will also be very interesting to compute the spectral function for the spinless fermion 
model. Considering several excited states and inserting the corresponding Bethe roots into 
the determinant representations for the form factors (3.12) and the scalar product (3.5), one 
can explicitly compute the spectral function for a finite system. The method developed in 
the study of the dynamical structure factors for the XXZ chain (see [21-26] for example) 
will also be useful in computations of the spectral function of the spinless fermion model. 
Comparison with the arguments in [27] is also interesting. 
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